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We study the various quantum aspects of the N = 2 supersymmetric Maxwell 
O 1 Chern-Simons vortex systems. The fermion zero modes around the vortices 

will give rise the degenerate states of vortices. We analyze the angular mo- 
■ mentum of these zero modes and apply the result to get the supermultiplet 

5— i 

structures of the vortex. The leading quantum correction to the mass of the 

vortex coming from the mode fluctuations is also calculated using various 

methods depending on the value of the coefficient of the Chern-Simons term 

k to be zero, infinite and finite, separately. The mass correction is shown to 

vanish for all cases. Fermion numbers of vortices are also discussed. 
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I. INTRODUCTION 



Abelian gauge theory in 2+1 dimensions with the Chern-Simons (CS) term has 
attracted much interest. Matter fields coupled with this term is believed to describe the 
anyons with the fractional spin and fractional statistics. Such important planar phenomena 
as the high T c superconductivity and the fractional quantum Hall effect has added more 
interest in the field theory models with CS term. The characters of allowed solitons are 
also affected by the presence of the CS term. As is well known, the usual 2+1 dimensional 
Abelian Higgs model supports only electrically neutral vortices as topologically stable soliton 
solutions ||. On the other hand, the CS term makes the vortices @] electrically charged, 
which are (extended) anyons J5J. We have quite rich vortex structures depending on whether 
the matter fields are relativistic || or non-relativistic j7| and whether we have more than 
one CS fields ||. In this work, we are mainly interested in the case with the relativistic 
matter fields coupled with the gauge field with both Maxwell and CS term in general. As 
a special limit of this general model, we get the Abelian Higgs model and the 'minimal' 
Chern-Simons Higgs model (i.e., without the Maxwell term in the action). 

With some special choice of the scalar potential in (2+l)-dimensional gauge models Ref. 
||, one can obtain interesting limiting theories in which the minimum energy static soliton 
solutions satisfy first-order differential equations, called the Bogomol'nyi [|TU[ or self-duality 
equations. This special potential becomes a specific scalar quartic potential for the Abelian 
Higgs model, while in the case of the 'minimal' Chern-Simons Higgs model a specific sixth- 
order potential form ||[ll|. The appearance of self-dual structures for certain special Higgs 
potentials can be ascribed to the extended supersymmetry Hl2-l^|. Requiring an N = 2 
supersymmetry guarantees this special form of the potential. There also exists an N = 1 
super symmetric model which produces exactly the same bosonic part of the Lagrangian as 
that of the N = 2 model. The fermion number is, however, not preserved in this case. 
We will mainly consider the model with more symmetry, i.e., the model with the N = 2 
supersymmetry. 
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A remarkable feature with these self-dual systems is the existence of static multi-vortex 
solutions which represent static configurations of vortices with unit flux without any inter- 
action energy between them. This interpretation is supported by counting independent zero 



modes |IT| , |i"5f to the boson fluctuation equations in the background field of a particular soli- 
ton solution. These bosonic zero modes are related with the collective modes of the solitons 
and play an important role in understanding the dynamics of the slowly moving vortices 
Of. 

Fermions around the vortices also have zero modes. For the models under study it is 
found that all fermion zero modes around the general multi- vortex background are closely 
related to the corresponding bosonic zero modes. The N = 2 supersymmetry is crucial |14| . 



They are very important in the quantum study of the models, representing the degeneracy 
of the soliton states (in contradistinction to bosonic zero modes which become collective 
coordinates) [|17]]. In supersymmetric models in particular, they account for the soliton 
supermultiplet structure [[Tj|. 

The organization of this paper is as follows. In Sec. II, we review the properties of the 
N = 2 supersymmetric Maxwell Chern-Simons system; The particle contents, the proper- 
ties of the vortices, the symmetries and the corresponding algebra of the Lagrangian are 
reviewed. In Sec. Ill, we will quantize the fluctuating modes around the self-dual vortices. 
By analysing the angular momentum fermion zero modes, we describe the multiplet struc- 
tures of the vortices. Spin contents of the degenerate supermultiplet of the vortex states 
are also calculated. Then, in Sec. IV, we calculate the mass correction to the vortices. We 
will do this for the value of k to be zero, infinite and finite, separately. Section V contains 
the summary of our work and discussions. Some technical details related with the spin 
assignment, supermultiplets and phase shift analysis are described in the appendices. 
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II. SUPERSYMMETRIC MAXWELL CHERN-SIMONS THEORY 



The Lagrangian for the Maxwell-Chern-Simons system with N = 2 supersymmetry is 



given by fH| 



C = C B + d 



(2.1) 



where 



and 



C 



B 



-\F, U F^ + -^ X F, U A X - \D^\ 2 - l -{d,Nf 
1 



(e|0| 2 + kN - ev 2 ) 2 - e 2 N 2 \<j)\ 2 , 



(2.2) 



(2.3) 



Here, = <9 M — ieA^ is the covariant derivative, N a real scalar, <ft a complex charged 
scalar, and ip (%) is a complex charged (neutral) 2-component spinor. Our metric tensor r] fiu 
has the signature (— , +, +). We will choose the 7-matrices as 7 M = (a 3 , ia 2 , 

When the coupling strength k for the Chern-Simons term becomes zero, the above La- 



grangian reduces to the N = 2 supersymmetric abelian Higgs model |T2[ . The scalar poten- 
tial in this limit allows only the symmetry broken vacuum. In another extreme limit of very 
large k (with the ratio — fixed), the neutral scalar field N (spinor field x) can be represented 
in terms of the complex scalar field (spinor field ip) as 



N = e 

K 



2 2\ 
— V ) 



X = —>/2e(f>*^, (2.4) 



and the Lagrangian becomes the supersymmetric extension of the minimal self-dual Chern- 



Simons Higgs model given in Ref. [13 

,(2) K 



OS 



,2\2 



+i^D^ - -(3|0| 2 - v 2 )^. 



(2.5) 



The above theory posesses various kinds of symmetries and their corresponding currents. 
The energy-momentum vectors related with the translational symmetry is given by 

P» = J d 2 xQ°» (2.6) 

with the energy momentum tensor given by 

e M „ = Ffj,\F u x + D^*D v (j> + D u <j)*D^ + d^Nd u N - i^D^ - ixi{^ u)X + V^C. (2.7) 

This energy momentum tensor is gauge invariant corresponding to the translation supple- 
mented with an appropriate gauge transformation. The generators for the Lorentz symmetry 
can be also found. For example, the canonical angular momentum operator J is given by 

J = J (fxfajxtPj - -(M + xx)\ 
= J b + Jf (2.8) 

with the contribution from the bosonic fields 

J B = J #xe ij x i [Ef i Nd i N + D°<p*DJ<f) + D J <P*D°<P + F ok F j k ], (2.9) 
and that from fermions 

J F = -iJ (fxe^x^D^ + xffyx] ~\j d 2 x{^ + %x)- (2.10) 
The theory in Eq.( [2.1|) possesses the following supersymmetry: 

S V A„ = HfffnX - XTnV), 

5, ? = v / 2#, 5 V N = i(XV ~ VX), 

5 11 X = Yv(d»N-f,)+ir ] G, (2.11) 

where 

F = eN<j), f ll = ~e lwX F vX , G = e|0| 2 + kN - ev 2 . (2.12) 



Here the spinor parameter r] should be taken as being complex Grassmannian. The corre- 
sponding supercharges are 

Q = V2j d 2 x[(D^)*Yl°^ ~ v / 2iF* 7 V 

^{d.N + f^l'x-G^x}. (2.13) 

The algebra of these with the supercharges are given by 

{Q a , Q?} = TrtfPp - 25?ev 2 <S>, (2.14) 

with $ = / d 2 xB. 

The particle spectrum will form representations of the above symmetry algebra. General 
structure of supermultiplet and the spin assignment will be described in the appendix. In 
the case of k = where there is only symmetry broken vacuum, all the particles are massive 
with the same mass v^2|ei>|. We have two (massive) vector modes with spin 1 and —1, two 
real scalar modes of spin and four spinor modes with two of spin 1/2 and the other two of 
spin —1/2. For the fermions, the sign of the mass term will depend on the spin. Two sets 
of the N = 2 supermultiplets are formed. One set is with one of spin 1, two of spin 1/2 and 
one of spin 0. The other is with one of spin —1, two of spin —1/2 and one of spin 0. 

For the Chern-Simons theory, the potential allows both symmetry unbroken and symme- 
try broken vacua. In the symmetry broken vacuum, we have four degrees of freedom with 
the equal masses 2e 2 v 2 /\k\ forming one set of N = 2 supermultiplet. The spin contents for 
k > are one of spin —1, two of spin —1/2 and one of spin 0. With k < 0, all the spins in 
the supermultiplet will change signs. In the unbroken vacuum sector, all the four modes are 
massive again with masses equal to e 2 v 2 /\k\. These are split into two supermultiplets. One 
supermultiplet consists with each of spin and spin 1/2, and the other with each of spin 
1/2 and I. If the sign of k becomes negative, the signs of all the spins are also changed. 

In the general case with finite value of k, there are also two degenerate ground states 
i.e., a symmetric one where <p = 0, N = ev 2 / k and an asymmetric one where \(f>\ = v, N = 0. 
The particle contents in the symmetry unbroken phase are the complex scalar and the 
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Dirac fermion ip with the mass e 2 v 2 and the neutral scalar N, the gauge field and 
the Dirac fermion \ with another mass equal to They form four N = 2 supermultiplets 
with two degrees of freedom. One supermultiplet consists of spin 1 and 1/2 and the other 
three consist of spin 1/2 and 0. In the broken phase, we still have two mass scales |(k 2 + 
4e 2 t> 2 ± \J ' k 2 (k 2 + 8e 2 v 2 )). The fields corresponding to these mass eigenstates are obtained 
as some combination of the original fields. We have two N = 2 supermultiplets. The spin 
contents in the supermultiplet with the mass m 2 + = |(k 2 + 4e 2 t> 2 + ^Jk 2 (k 2 + 8e 2 v 2 )) are one 
with spin 1, two of spin 1/2 and one of spin 0. The spins in the other supermultiplet with 
m 2 - = \{.k 2 + 4e 2 f 2 — J k, 2 (k 2 + 8e 2 v 2 )) will be that of the m + with all the signs changed. 
The brief analysis of the above supermultiplet structures is done in the appendix. 

Now, let us briefly review the structure of self-dual vortices in the Maxwell-Chern-Simons 
theory. The bosonic Lagrangian Lb m Eq.( [2.5|) is enough for the classical solution. In 
this theory, there are two degenerate ground states as mentioned above. It is known that 
topological solitons exist in the asymmetric phase with the asymptotic behavior 

N(r)->0, \<f>(r)\->v as r -> oo (2.15) 

and a quantized flux $ = ±^n (n: positive integer). Nontopological solitons also exist in 
the symmetric phase. We will consider only the topological vortex for simplicity. These 
vortices satisfy the field equations. Especially, they satisfy the Gauss law constraint 

diF i0 + kF 12 + e J° = (2.16) 

with J° = -i((f)*D°(j) - D°(j)*4>). Integrating over the whole space then tells us that a 
configuration with the magnetic flux $ carries the electric charge Qe = / d 2 xJ° = — -<&. In 
this theory it has also been shown || that the energy of the configuration is bounded from 
below by the relation 

E > ev 2 \$\ = 27rv 2 n, (2.17) 

and is saturated if the configurations satisfy the following 'self-duality' equations 
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{D 1 ±iD 2 )<j> = D ± <f> = 0, 



F 12 ± (e|0| 2 + kN - ev 2 ) = 0, 



A°^N = 0, (2.18) 



together with the Gauss law ( 2.16 ). The upper (lower) sign corresponds to a positive (neg- 
ative) value of the magnetic flux $. Whenever we need the explicit choice of the self-dual 
background field configuration, we will choose the upper sign corresponding to the vortex. 

In the case of k = 0, we may consistently set A = N = and Eq.( |2.1£| ) will become the 
self-duality equations for Landau- Ginzburg vortices [ ID | . On the other hand, in the case of 
k — > oo, we have instead 

A ' = &W (2 - 19) 

and the self-duality equations reduce to those of Ref. ||, viz., 

D ± <P = 0, 

F 12 ±^-\<p\ 2 (\<p\ 2 -v 2 ) = 0. (2.20) 

We can take the spherical ansatz for those classical vortices with vorticity n on top each 
other 

= f(r)e me , eA l = e^(a(r) - n). (2.21) 
The functions are related as 

a(r) =r— In f(r), (2.22) 
dr 

and can be solved using the Bogomolyi equation. Now consider the angular momentum in 
the presence of the classical vortices given by the spherical ansatz. With the Gauss law now 
containing the fermion charge density in J°, the angular momentum from the bosonic fields 
in (|2.9|) can be written as 

J B = J dPxEijX* [d°Nd j N + + &(j>*D°(j> + F ok F j k 

+A j (-d k E k + KF 12 + e^°i})}. (2.23) 
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The leading contribution to the angular momentum in (|2.23|) for the vortices will be from 
the background classical vortex configuration J c i and from the fermion zero modes AJ - 
Jb ~ Jci + AJ . For vortices with the spherical ansatz, these are given by M 



Jd = y-n 



(2.24) 



and 



AJ = - / d 2 xa(r)^^). 



(2.25) 



The fermionic field contribution to the angular momentum in Eq. fl2.10p under the vortex 
background becomes 



J F = I d 2 x[i)\-id e + a(r) - -a 3 - n)i/j + x\~id e - -x^xl- 



[2.26) 



The total contribution from the fermion modes will then become 



J F + AJ = J d 2 x{ip\ X ] ) 



-ida 



f a 3 o\ ( n 0^ 

a 3 



V 00 / 



[2.27) 



This will be used in the next section. 



III. QUANTIZATION 

Now, we want to quantize the theory in the soliton sector. The general procedure of the 
quantization of fields around the soliton has been well developed We decompose the 
fields around the classical vortex configuration as 

$ = $ d + 5$. (3.1) 

The field 5$ is the fluctuating modes around the classical vortex configuration $ c j. Plugging 



eq.( |3.1| ) into the Lagrangian, we have the Lagrangian in the form of 

£ = C c i + C(2) + Cinf. 
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Here, L c \ is the same as the bosonic Lagrangian in Eq.( [2.2|) except that all the fields are the 
classical vortex configuration. This is the tree level contribution to the Lagrangian coming 

(2) (2) 

from the classical vortex configuration. £( 2 ) — £>b + * s the quadratic piece in terms of 
the fluctuation fields. For general value of k, the quadratic part of the bosonic fluctuations 
is given by 

Cf = --SF^JF^ + - a k^ x 5F^5A x - \D^<P\ 2 - e 2 <L4j|0| 2 

-ieSA^D^ + (j>*D^8(j>) - \{d^5N) 2 - ^(0*50 + 050* + ^6N) 2 

-e 2 (|0| 2 + -N- v 2 )\5(j)\ 2 - e 2 (^) 2 |0| 2 - e 2 AT 2 |50| 2 - e 2 AWV(0*50 + 050*), (3.2) 
e 

and that of the fermionic fluctuations by 

CP = ifyfD^ + ixi^d^x + 

-iV2e{tpx<P - X*P<P*) + eN^ip. (3.3) 

All the bosonic fields above is the classical background. Terms in Li nt are the higher order 
interaction terms. 

The equations of motion for the fluctuating fields around the self-dual vortices can be 
obtained by varying the quadratic piece of the Lagrangian in the above. First, the equations 
of motion for the bosonic fluctuations are 

(d^ - 2e 2 \(f)\ 2 )5A u - d^d^A" - ne^d^Ax - ie{<j>*D v '50 + 50*5» = 0, 
(-d 2 + D_D + - 2ieA°d )6(f) - ie(D^<f))5A + - Ke<j)5N + 2eA°0(5A° - 6N) 

-ie^d^A" - e 2 (j)((j)*6(j) + 050*) = 0, (3.4) 
(d^ - 2e 2 |0| 2 - k 2 )5N - e(n + 2eiV)(0*50 + 050*) = 0. 

Time independent solutions of these equations are the zero modes. The bosonic zero 
mode fluctuations may also be obtained by considering the variation of the self-duality 
equations ( |2.18| ) around the given classical vortex configuration as given in Ref. ||15|| . Among 
the zero modes, we eliminate those related with the gauge transformation by imposing the 
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gauge fixing condition. The index theorem or its variants can be used to count the number 
of bosonic zero modes satisfying the equations of motion and the gauge fixing condition. The 



bosonic zero modes for the k = pOfl , those for the k = oo case |Tl| and those for a general 



value of k were studied ||15|| . The results are that in the background of a topological vortex 
configuration with vorticity n, there exist 2n bosonic zero modes for any value of k. They 
correspond to the collective coordinates associated with the vortices. The quantization of 
the bosonic zero modes will give rise to the excitation of the collective coordinates e.g., the 
momentum. This is consistent with the interpretation of these zero modes as being related 
to translation of individual vortices. 

The general time dependent modes can be quantized as usual to describe the bosonic 
particle excitations around the vortices. 

We now turn to the fermion fields. The Dirac equation for fermion fields around the 
vortex is 



idn 







(A 


— Hp 









(3.5) 



where the Hamiltonian Hp is given by 

/ 

H F = 



-i7°7 ■ D + e(-j°N + A 



V 



—i 







v / 2ie7°0 
— ry°7 • V — 7 k 
\ 



V 



o 

V2e^ 

d + 
in 



(3.6) 



D + 

D- 2ieN 
V2e(f)* -m 
-V2e<p* 9_ 

The background fields in the above are for vortex configuration corresponding to the upper 
sign in Eq.( [2.18 ). 

Solutions of the Dirac equation of fermions around the vortex in the above equation can 
be decomposed as 



(3.7) 
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Here are the zero modes and £) w b^ +UJ (E w ($*&-u) are the positive (negative) energy 



solutions. The fermionic zero modes are analyzed in Ref. |TJ| using the index theorem and 
also the relation between the bosonic zero modes and the fermion zero modes. There are 
2n fermion zero modes around the vortex configuration with winding number n as in the 
case of bosonic zero modes. This is deeply related with the N = 2 supersymmetry of the 
theory fll4 |. 

The quantization of the zero modes of the fermions will be relevant to the multiplet 
contents of the vortices. To do this, we need to know the angular momentum of the zero 
modes. The quantum mechanical angular momentum operator J can be read from the field 
theoretical expression in Eq. ([2.27|) as 



J 



V 



-»J u 3 2 " 




/ 



(3.8) 



—ieijdj — | CT3 — n 

—ieijdj - |cr 3 

It is straigtforward to show that the angular momentum operator J commutes with the 
Hamiltonian Hp. We decompose the modes into angular momentum eigenstates. The 
general mode with the angular momentum quantum number j is given by 



\XlJ 

The Dirac equation for this mode becomes 
/ 



^ h 1 (r)e i ^ +n ^ e 
h 2 {r)e i{j -^ e 



\ 



(3.9) 



h x {r) 
h 2 (r) 

h 4 (r) 



d r + 





q+j+l/2 



r 



-V2evf 




V 



V2evf 




2ieN 

—in d r 

-V2evf d r + i^ 





V2evf 



./-1/2 



IK 



fti(r) 



h 2 (r) 
h(r) 



J 



(3.10) 



We first consider the fermion zero modes in the k — > limit. In that case we have two 
sets of equations. 
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-V2evf d r + ^ 



( h 2 (r) ^ 



(3.11) 



and 



( dr+ a ±mI1 \ ( h J^\ 



\ 



\/2evf d T 



j- 1/2 



(3.12) 



The second set of equations ( |3.12j ) does not allow any normalizable solution. The first set 



of equations (|3.11|) can be combined as 



r H 



h 9 



V2evf 



0. 



(3.13) 



The asymptotic behavior of h 2 {r) in ( |3.13|) will be proportional to e ^ evr . Near the origin, 
h 2 behaves as 



h 2 ~ r n {A x r j ~^ + A 2 r~ {j ~^). 



(3.14) 



We expect single solution matching the boundary condition at infinity with two free param- 
eters at the origin. To get the regular solution at the origin with two free parameters, j is 
restricted to half integral values in — n + \ < j < n — |. The value j = n + | is discarded 
since corresponding ip^ 



Vi ~ (A ir 2n + A 2 )t 



(3.15) 



has bad behavior at r = 0. For each of the above 2n solutions of ipi (h 2 ), the function if)^ 



(hs(r)) is determined through Eq. (|3.11|) . Hence we have 2n independent zero modes. This 



result agrees with that in Ref. [14| based on the index theorem. 

For the general value of k, we can get two coupled 2nd order differential equations of h 2 
and hi from (|3.10| ). 



d 2 r + 



d r (j ~ 1/2) S 



2e 2 v 2 f 



ho 



V2evf 



+ IkHa = 0, 



(3.16) 



and 
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d 2 + dr_ (j 1/2) 2e 2 v 2 f 2\ ^ _ - (/{ + 2ejV)2e 2 l; 2 / 2 ( ^2_ ) = Q ^ 



V2evf 

The remaining functions /i3 and h\ can be determined by h 2 and ^14 through the 3rd line 
and 1st line in the Dirac equation in (|3.10| ). The general solutions of Eqs.( |3.16|) and (|3.17|) 



are expected to have four free parameters, to be adjusted to the boundary conditions at 
the origin and infinity. Near r ~ 0, the regularity of the solution gives us only three free 
parameters for — n + | < j < n — |. The leading orders in the power series expansions in 
that range of the angular momentum are given by 

h 2 ~ r n {A ir j -^ + A 2 r~ (j -^), h 4 ~ Br^l (3.18) 

If the angular momentum is out of the above range, then we have at most two free param- 
eters. In the asymptotic region, among the four parameters, two will conrrespond to the 
unphysical divergent solutions and only two free parameters will show up in the convergent 
solutions as 

h 2 ~ C ie - mir + C 2 e- m2r , -K ~ mi2 ~ 2eV C!e- mir + m 2 2 -2e 2 v 2 c ^_ m2r ^ 

in k k 



where m 2 and m 2 2 are eigenvalues of the following mass matrix 

(3.20) 

Matching the solutions to have the regularity in (|3.18|) at the origin and the integrability in 



( 2e 2 v 2 k 2 ^ 
yn 2 + 2e 2 v 2 2e 2 v 2 j 



( 3.19 ) will then leave us only one free parameter that comes from the homogeneity of the 
differential equations. In other words, we have 2n solutions, one for each j in the range of 
—n + \ < j < n — |. Note that we do not expect any solution if j is not in the above range, 
since we have less parameters in the power series solutions near the origin. Specifically, for 
n — 1, we have two modes with j = ±|. 

Based on this analysis, we quantize the theory. For simplicity, we consider the case of 
single vortex with the winding number n = 1. Multivortex case can be similary done when 



they are widely separated. The quantization b and d in Eq.( [3.7|) (with their conjugates) for 
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nonzero modes are the same as that in the vacuum sector and these modes describe the 
fermions around the soliton. We have two fermion zero modes with the angular momentum 
±|. The quantization for these modes will be 

{a i ,a}} = 6 ij , (i,j = l,2). (3.21) 

The subscript 1 represents for j '< = — | mode and 2 for j = |. According to Jackiw and 
Rebbi's interpretation ||17| , the soliton states will be degenerate due to the fermion zero 



modes and the quantum multiplet structure and the spin contents of the vortices will form a 
representation of the algebra relations in equation Q3.21J) . We will then have four degenerate 
soliton states from the fermion zero mode algebra in (|3.21|) . 



| >, |+— >=0i| >, |— + >= <4| >, \ + + >=a\a\\ >. (3.22) 

The algebra in Eq. (|3.21|) is not the same as the N = 2 SUSY algebra with the central 
charge described in the appendix. The discrepancy comes from the fact that only the mode 
corresponding to j = — | can be obtained by supertranslation of the vortex configuration. 
To see this, note that we can always get one fermion zero modes by the supersymmetry 
tranformation in Eq. ( |2.11| ) to the classical bosonic vortex background. We then get one 
fermion zero modes ^ proportional to 

ipl = 2V2F, ^ = V2iD_(p xt = 2 ^G, xi = ~2d^N (3.23) 

It is straightforward to check that the spin of this zero mode is j = — | using the angular 
momentum operator in (|3.8| ). Hence this is the zero mode corresponding to a\ in (|3.21| ). We 



now write the supercharge in (|2.13|) in two component form 

f V2(D 4>* - iF*)^ - i(d N + f - i G) X] + yfiiP-W^ - i{d + N + f + ) Xl ^ 
v V2(D 0( p* + iF*)ij i - i(d Q N + f + iG) Xi ~ V2(D + <f>Y^ + z(«9_iV + /_) Xt ) 

(3.24) 



Q = Jd 2 x 



Note that the down component of the supercharge with the background self-dual bosonic 
fields vanishes. The nonvanishing upper component of the supercharge becomes 
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2 J d 2 x (-y/2iF*^ - G XT + ^(D^)*^ - if +X )j • (3-25) 

The algebra of the supercharge will then be that of N = 2 with the central charge in the 
appendix. The nonvanishing supercharge is in the form of = J ^x^f^. From the 
orthogonality of the modes, this is proportional to a x and so the operator a 2 corresponding 
to the other zero mode anticommutes with the supercharge. In other words, among the two 
independent supertranslations to the self-dual background configurations, one from Q± acts 
trivially and only the other one from will give us the fermionic zero mode corresponding 
to j = —\. The other zero mode is not obtained from supersymmetry. The algebra between 
d\ and a\ is the same as the N = 2 SUSY algebra with the central term in the appendix 
and will be realized as a doublet state. On the other hand, the doublet representation of the 
algebra from a 2 and a\ will transform as a singlet under the N = 2 SUSY. Hence the above 
four degenerate solition states will form two sets of N = 2 supermultiplets rather than one. 

One supermultiplet will be by | > and | H — > with angular momentum J c i and J c i — 1/2 

respectively. The other one is by | — h > and | + + > with angular momentum J d + 1/2 
and J c i respectively. Here J c \ is the leading contribution to the angular momentum from the 
classical bosonic field configuration of self-dual vortices. 

We now calculate the fermion number of the soliton by taking the expectation value of 
the fermion number operator for the degenerate soliton states. 

\\ J d 2 xMx)\ ¥(*)]! ± ± >= ±1 ± \ + V (H F ) (3.26) 

The constant pieces differing on the vortex structures are from the fermion zero modes and 
t)(Hf), the so-called 77-invariant, given by 

V(H F ) = W-E)- (3-27) 

is from the nonzero modes. 



< ±± 



For the Landau- Ginzburg model, i.e., k = case, there exists a constant matrix 
that anticommutes with the Hamiltonian. This matrix matches the positive 



( 7° ^ 



- 7 U 
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energy solution with the negative energy solution with the norm and the density of states 
preserved, making the value rj to be 0. This can be also shown easily by direct evaluation 
of rj using the method in Ref. for example. Then from Eq.( |3.26| ), the four degenerate 
states of the vortex (|3.22| ) carry fermion number —1, 0, 0, 1, respectively. 

The Hamiltonian (|3.6|) for the general value of k no longer has such a structure. Niemi and 



Semenoff developed a method to calculate rj for such a general case |2T| which are based on 
the works of Atiyah, Patodi and Singer . Following them, let us introduce one parameter 
family of Hamiltonian H(t) which interpolates the Hamiltonian Hq = Hf(k = 0) (when 
t = — oo) and the Hamiltonian Hp (when r = oo). The Dirac operator D T = 27 (d T — H T ) 
is defined on the extended manifold A4 of R 1 x D 2 where R 1 = {r} and D 2 is a disk of 



radius R in the usual x-y plane. The result is [pi]] , in the R — > 00 limit, 

- l V {H F ) = Index(P T ) - ~ V (H ) + ^(Re(P)) (3.28) 

Here, Index(D T ) is the index of D T in the extended manifold Ai and P is the operator 
defined by projecting the operator D T onto the boundary of the disk D 2 for each value of r. 
For our case, it becomes 



ReP = -i-f°d T + i 



^ o\ cos 6 + o"2 sin 9 ^ 



d e (3.29) 



0i cos 9 + o" 2 sin 9 

which is almost free equation. One can easily see that 7/(ReP) = and i](Hq) = 0. Hence 
we get tj(Hf) = — 2Index(P T ). The index of D T is hard to evaluate. Since the index is an 
integer, r] is an even integer in general. But we expect its value to be zero, since 77 is shown 
to be zero in the above when k = 0. In other words, there is no contribution to the fermion 
number from the nonzero modes. Then the fermion number of the vortex is the same as 
that in the k = case. 



IV. QUANTUM CORRECTION TO THE MASS OF THE VORTEX 



The quantization of the nonzero modes in Eq. ( |3.7| ) will correspond to the particle creation 
and annihilation around vortices. The leading mass correction comes from the quantum 
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fluctuation of these modes. We will get this correction by comparing the bosonic and 
fermionic modes. 

The quantization of the nonzero modes in Eq.( |3.7| ) will correspond to the mass correction 
of the vortex and particle creation and annihilation around vortices. The mass correction 
comes from the sum of those mode contributions |^3[^3|. We will get this correction by 
comparing the bosonic and fermionic modes. First, let us calculate the quantum correction 
to the mass of the vortex. The leading quantum correction to the vortex mass comes from 
the quantum fluctuation modes given by 

,dns{X) dnp(X) . 



AM = Y^^b - = / d\\f\{- 



-)■ 



(4.1) 



dX dX 

The frequency ojb {^f) are eigenvalues and n B {n F ) are the number of states upto eigenvalue 
A of the bosonic (fermionic) fluctuating modes. 

We first evaluate the mass correction for the usual Maxwell theory (k = 0). The equations 
for the fermionic modes around the self-dual vortex with the positive frequency up becomes 



110 p 







< D F \ 


(u\ 






\ -D\ J 





(4.2) 



where 



u = 




v = 


M 




v Xl J 







and the Dirac-like operator Dp is defined as 

( 



\ 



-V2ec 



0_ 



\ 



(4.3) 



Here the bosonic fields and Aj in the covariant derivatives are the classical background 
fields of the self-dual vortex. Apply the Dirac like operator Dp to the equations of fermions 
in Eq.Q to get 



' D F D F 



D F D F 



(4.4) 



The contribution from the fermionic modes to the soliton mass is given by 

Y^v F = J2 u u + J2 u v- (4.5) 

Note that if (U, V) is the mode corresponding to uj f then (U, —V) is the mode corresponding 
to —up. Let us represent -jf (-yr-) as the density of states of the operator of D F D F (D F D F ). 
Then the density of the states for the fermion modes in Eq. (f4.2|) is half of those of the 2nd 
order equation in Eq.( [4.4| ) 

dn F {\) 1 dn+(X) dn-(X) 

-^jT = 2 ( ^jr + ^)T ) > (4 - 6) 

since only half of the solution of the above equation correspond to the positive frequency 



solution in Eq.(4~2). 

We now turn to the bosonic fluctuations. The equations of motion in Eq. ( |3.4[ ) for k = 
become 

(d^ - 2e 2 \(j)\ 2 )5A + - d + {d^A* + ie(<j)*8<j) - <p5(j)*)} + 2ie(£L0)*<ty = 0, 
(-d 2 t + D_D + - 2e 2 \(f)\' 2 )5<f) - iecf) {d^SA* + ie(<j)*5<j) - <j)8<j)*)} - ieD^(j)5A + = 0, 

(V 2 - 2e 2 \(j)\ 2 )5A + j t [WidA' + ie(<p*5<p - <p5<p*)} = 0, 

{d^ -2e 2 \(j)\ 2 )5N = Q. (4.7) 

The gauge field is written in terms of 5A + = 5Ai + i5A 2 for convenience. The classical 
background fields and A4 appearing in the above equations are the configuration for the 
self-dual vortices satisfying the self-dual equations (self-dual) with k = 0. Among the 
fluctuations satisfying the above equation (|4.12|) , we have to subtract those fluctuation 



modes corresponding to the gauge transformation. For this purpose, we choose the physical 
gauge condition as follows. 

VM* + te((f)*5(f) - <f)6(f)*) = 0. (4.8) 

The above Eq. ([4.7|) for SA° with this gauge condition then requires that 

5 A = 0, (4.9) 
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since the operator V 2 — 2e 2 |0| 2 is negative definite. Other equations for the bosonic fluctu- 
ations can then be written as 



(d^ - 2e 2 \<p\ 2 )5A + + 2ie(D_0)*<ty = 0, 
-d 2 + D-D+ - 2e 2 |0| 2 )50 - ieD^(f)5A + = 0, 
(d^ -2e 2 \<j)\ 2 )5N = 0. 



(4.10) 
(4.11) 
(4.12) 



To find the equations for the fluctuation modes corresponding to the gauge transformations, 
we take an infinitesimal gauge transformation 5q5A + = <9 + A and 5gS4> = ieA(p to equations 
( |4.12j ). We then get the equations for the gauge transformation modes satisfying 



-2e 2 \cj)\ 2 )A = 0. 



(4.13) 



We have to subtract the degree of freedom satisfying this equation. Since this equation is the 
same as that of the real scalar field of 5N in Eq. ( |4.12| ), the contribution of the fluctuation 
corresponding to the gauge degree of freedom A satisfying Eq. ( |4.13| ) and that of the neutral 
scalar 5N in Eq.( [4.12D cancels out in the contribution to the correction of the mass. Hence 
the only bosonic contribution is from the complex fields 5A + and 5<f) and the vortex mass 
correction from bosons becomes 



To get the density of the states we write the equations for 5A + and 5(f) in Eqs. 
matrix form as follows 



(4.14) 



in the 



d 



( *\ ( 



D_D + -2e 2 \(j)\ 2 \/2ei 

V2e<J)* di 2 - 2e 2 \<f)\ 2 
( \ 



V 7§ M + J 



(4.15) 



(4.16) 



where the Dirac operator Dp is defined in Eq.( |4.3j ). From this we see that the density of 
states for the bosons are n + . With the result for fermions in Eq.(J4.6|), this gives that 
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driB^X) dn F (X) 1 ( dn + {X) dri-(X)^ , , 

, YiL.lt) 



d\ dX 2 y dX dX 

and hence the mass correction becomes 



AM = J2 - Yl ( 4 - 18 ) 



Udx( d ^- d -^)VX. (4.19) 
2 J y uA clA / 

We can evaluate the integration by calculating the density of states from the phase shift. 
This method is summarized in the appendix C. and we get zero for the value of the above 
integration. 

Evaluation of the integrand can also be done with the help of the index formula [f20|1 . 

I{ , ) = Tt (' * ) = r«* *h£)Y (4.20) 

\z + D F D F z + DpD F ) Jo z + X\ dX dX J 

The index is easily calculated and the result is 2n. For the index to be independent of z, 
the integrand in the above should be 

dn + (X) dn_(X) 



dX dX 



2n5(X). (4.21) 



By plugging this result into the Eq.( [4.19 ) we see that the mass correction to the vortex in 



N = 2 model at one loop level vanishes. 

AM = 0. 



This result agrees with Ref. p5[ . 

As a check for the independence of the choice of the gauge, let us describe this in the 
covariant gauge. We choose the background gauge by adding the following gauge fixing term 
to the Lagrangian. 

C g .f = ~\ R<L4" + ie{<P*5<P - <M0*)} 2 • (4.22) 

The equations for the modes in this gauge are then obtained from the quadratic pieces of 
the Lagrangian by adding . 
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+ LP + Cg.f + Cgh- 



(4.23) 



The ghost Lagrangian from the gauge fixing is given by 



C gh = c(d^-2e*\<P\*)c 



(4.24) 



The equations for the fermionic modes are the same as before. For the equations of the 
bosonic fluctuations, those for 8A + , 8(f) and SN are the same as before in ( |4. 1 0|) and ( [4.111) . 
The equations for 5A° which is dynamical in this gauge and the ghost field c are exactly 
same as that of field 5N in ( [4.12| ). Hence the contribution from the bosonic fields and ghosts 
is given by 



(4.25) 
(4.26) 



We have used the fact that the contribution to the mass from the ghost fields and that from 
the field 5N and 5 A cancel out since the equation for the ghost fields are the same as those 
for fields 5N and SA° in Eq.( 4.12 ). We have shown that the only contribution comes from 
5A + and 5(f>. Since equations for these fields are the same as those in the Coulomb gauge, 
the remaining arguments are the same as before and so we get the same result for the mass 
correction for the vortex. 

Now consider another extreme limit of the Lagrangian with k — >• oo. The equation of 
motion of the fermion modes are 



D 



(oo) 



(4.27) 



with 



D 



( 2i e -\d>\ 2 



(oo) 



V 



(4.28) 



D_ -2^(2|0| 2 -t; 2 ) 
The equations of bosonic fluctuations can be obtained in a straighforward way from the 



Lagrangian ( |2.5|) . 
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Ke^dJAx - 2e 2 \<f)\ 2 5A» - 2e(0*Z>50 + <50*Z>0) = 0, 



(4.29) 



and 



Dl5(j> + ie(25A°D + d 6A°)(f) = D 2 5<f) - ie{5A i D i <j> + D i 5A i 



- ^(9|0| 4 - 8v 2 \<b\ 2 + - £-0 2 50*(6|0| 2 - 4v 2 ). (4.30) 



We choose the gauge for the spatial gauge fields as 



V % 5A l + 2i— (2|0| 2 - v 2 )((f)*5<j) 

K 



—((j)*d Q 8(j)-^8(j)\ 

K 



(4.31) 



The zeroth component of Eq. flPDp and the gauge fixing condition ( PI ) can be combined 
into the complex equation (and its complex conjugate) 



3 2 

d-5A+ + 4^(2|0| 2 - v 2 )(j)*5(t) - 2i—\(P\ 2 Q = 2-<f>*d l 



K 



K 



(4.32) 



where 



Q = 5A° + 



+ 



K 



(4.33) 



The spatial component of Eq. ([4.29 ) can be written as 



d SA + + 2i— \cf)\ 2 5A + - 2-(P*D + 5(p + d + Q = 0. 

K K 



(4.34) 



We will show that the gauge fixing condition and the above equations of motion gives the 
SA° so that Q in the above equations is zero, and 5 A is fixed as 



5A C 



K 



(4.35) 



To show this, we first rewrite the equations for the scalar field (|4.30|) using the gauge fixing 
condition and the equation ( |4.32|) as 



d 2 5(p = D_D^ 



4-^(2|0| 2 - v 2 ) 2 54> - ieSA+(D_ 



+ie— (|0| 2 - v 2 )d_5A + - ie4>d Q. 



(4.36) 
(4.37) 
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Comparing this equation with the equation obtained by taking the time derivative of the 
equation ( |4.32| ) we get 



V 2 Q = V'(cL4 u + -U*5(p + 8<j>*(t>)) = 0, 

K 



(4.38) 



hence Q = 0. The equations ( 4.32p and (4.34) then become 



*d 5<f> = ^-d^5A + + 4i-(2|0| 2 - v 2 )0*<50, 



(4.39) 



and 



d n 5A A 



-2i—\(f>\ 2 5A + -2-0*D^ 

K K 



(4.40) 



And they can replace the set of Eqs. ( 4.30|) , ( |4.32| ) and ( f4.34|) . The imaginary part of 



the equation (|4. 39|) is nothing but the gauge fixing condition in Eq.( |4.31|) . We can also 
easily see that the two equations (|4.39|) and ([4.40|) with 5A given as in Eq.( f4.35|) imply 



the equation of motion of the scalar fluctuation in Eq.( [4.30| ). Hence the physically relevant 
bosonic fluctuation modes are described by the above two equations. 

To compare these equations with the fermionic modes we write the equations fl4.39p 
andflggg). 





1 


k SA+ \ 




( k 5A+ \ 


do 




2e 4>* 


= -D(oo) 


2e <j>* 




V 


H ) 







(4.41) 



Note that this equation has precisely the same form as that for the fermionic modes in 
equation ( |4.27|) . This means that the density of the modes as well as the spectrum of the 
bosons are equal to that of fermions. Hence the quantum correction for the mass of the 
vortex vanishes identically. 

Finally, let us consider the general case with finite k. The equations for the the bosonic 
fluctuation fields are given in Eq.( |3.4| ). We have to fix the gauge to eliminate those fluctua- 
tions corresponding to the gauge transformation. We choose the following background type 
gauge condition. 
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did A 1 + k8T + ie( 



-d n 5N, 



(4.42) 



where 8T satisfies the following equation. 



doSF - eijdiSAj - kSN - e(0*50 + 



0. 



(4.43) 



The equations of motion in Eq. (|3.4j ) with the help of the gauge fixing condition in equation 
( [4.43 ) can be written as 



(8 2 + D_D + - 2e 2 |0| 2 )50 - 2ieA°d 5(j) - ieD_<j)8A + - neSNcf) 

+2eA°(j)(5A - 5N) + ietufiT = 0, 
(d 2 - 2e 2 \(j)\ 2 )5A - (k 2 + d 2 )8N - e(« + 2eiV)(0*<50 + 50*0) = 0, 
(d M 0" - 2e 2 \(j)\ 2 )5A + + Kd + 8F + m{d + 8N + d 8A + ) + 2ie(LL0)*50 = 0, 
(8^ -k 2 - 2e 2 \(j)\ 2 )8N - e(n + 2eiV)(0*50 + 050*) = 



(4.44) 
(4.45) 
(4.46) 
(4.47) 



We compare Eqs. fl4.45] ) and (|4.47[) to get 



(d 2 -2e 2 \4>\ 2 )(5A° -5N) = 0. 



(4.48) 



Then 



8A° = 5N. 



(4.49) 



This relation reduces (|4.44 ) and ( |4.46| ) in the following form. 



-dl + D_D + - 2e 2 |0| 2 )50 - 2ieNd 5(j) + ien${8T + i8N) - ieD_ 
(-8 2 + V 2 - 2e 2 |0| 2 )A+ + iKd 8A + + k8+(8F + t8N) + 2ie(D. 



0. 
0. 



(4.50) 
(4.51) 



We have a set of equations (|4.50|) , (|4.51 ) and (|4.47|) which are supplemented by the gauge 
condition flP2D and (fl43|) or 



(8b + iK,)(6F + i6N) + id^5A + - 2e0*50 = 



(4.52) 



We now turn to the Dirac equation given in Eq. ( |3.5p . We can remove from this 
equation by using the first line of that equation 
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- doV T = D+tf>i - V2e<f)xp (4.53) 

Then we have 

+ D_L> + - 2e 2 |0| 2 )^ - 2ieNiP l - v^e(D_0)x T + V^ie/c^ = 0, (4.54) 

and 

(-«9 2 + V 2 - 2e 2 |0| 2 ) XT + «^oXT - V2e(D^)*^ l + md +Xi = 0, (4.55) 
together with the last line of the Dirac equation. 

d oXl - >/2e^Vi + d-Xi + iKXi = °- ( 4 - 56 ) 



Comparison of bosonic equations (|4.50|) , (|4.51|) and (|4.52|) with fermionic equations Q4.54|) , 
( |4.55| ) and (|4.56|) gives us the relation 



^ = S</>, XT = 4= <*A+, xi = -j=(8F+i8N). (4.57) 

Note that this identification was pointed in Ref. |l4j] for the zeromode case. Here we show 
that the identification holds also for nonzero mode. The number of fermion states np is 
simply related with the sum of the phase shifts of ip^, ipi, \h an< ^ Xi- m the Appendix C, 
we will show that the sum of phase shift of ij}^ and xi is equal to the sum of phase shift of 
ipl and Xv Therefore 

n F = + n Xv (4.58) 

For the bosonic case, the physical degrees of freedoms are 5(f), 5N and 5A + where gauge 
degree of freedom is subtracted. However Eq. ( [4.49 ) says that d^SN = iuSN (assuming the 
time dependence as e~ lut ) is nothing but 

ViSA* - inSF + ie((j)*5(j) - 050*) = -u6N, (4.59) 

which is gauge degree of freedom. So we assert that, as in the case of k — 0, we would 
subtract 5N instead of subtraction the gauge degree. Then bosonic degree of freedom is 
given by 5(f) and 5A + . Then 
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ns = nscp + usa 



■+ 



(4.60) 



and so we get 



Tip = TLB- 



(4.61) 



Therefore there is no mass correction 



AM = 0. 



(4.62) 



V. SUMMARY AND DISCUSSION 



We have studied various quantum aspects of the N = 2 supersymmetric Maxwell Chern- 
Simons theory. First we identified the mass spectrum, the spin contents and the supermul- 
tiplet structures of the particles both in the broken and the unbroken sectors. 

Then, we analyzed the vortex sector, the main subject of this paper. Starting from the 
canonical angular momentum we evaluated the leading quantum correction to the classical 
value of the angular momentum of the vortex coming from the fermion zero modes. For 
the supermultiplet structure of vortices, fermion zero modes play an important role. The 
algebra by the operators of the fermion zero modes around the winding number n = 1 
vortex is larger than that of the N = 2 SUSY algebra with the central charge. They provide 
two supermultiplets with the relative spin difference half, rather than single supermultiplet. 
This is in contrast with the case in the monopoles in the 3+1 dimension or kinks in the 1+1 
dimension. The fermion number of the vortex is also calculated. 

Leading quantum correction to the mass of the vortices is calculated separately depending 
whether k = 0, k = oo or finite k. The mass correction can be obtained either through the 
index theorem or by comparing the modes between bosons and fermions. In all cases, we do 
not see any mass correction. 

Self-duality is deeply related with the underlying supersymmetry. Here we have consid- 
ered only N = 2 supersymmetric model. For the model with N — 1 supersymmetry that 
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allows self-dual vortices is not treated. We expect mass correction in this model since we do 
not see any simple way of matching the bosonic and fermionic contributions. 

In the models with the self-duality, the mass of the vortex will be simply related with 
the magnetic flux of the objects at the tree level. An interesting question is whether this so- 
called Bogomolnyi bound will still be saturated at the quantum level. There are some models 



in 2 and 4 dimensions known to satisfy the Bogomolnyi bound at the quantum level fl26| , |18 



To show the saturation at the quantum level, Olive and Witten |26[ used the argument 
that the size of the supermultiplets for particles and solitons cannot change abruptly by 
perturbation. This argument seems not to be directly applied in our case. First, we don't 
have any self-duality for particles since we don't have any conserved charge for the particles 
in the broken sector. In the vortex sector, for the single vortex with the winding number 1 
to be specific, we have four degenerate vortex states from two fermion zero modes. These 
form two irreducible supermultiplets of size two with the Bogomolnyi bound saturated. If 
the bound is not saturated, we might have only one supermultiplet of size four. Whether the 
bound is saturated or not, the total size of the states remains the same. On the other hand, 
in those models in other dimensions mentioned above, the supermultiplets of the degenerate 
solitons form single irreducible representation of the superalgebra with the bound saturated. 
If the bound does not become saturated by the perturbation, we need more states, which 
is unlikely. This is the difference between our models and those in other dimensions. This 
difference is related to the fact that all the fermion zero modes in our case do not come from 
the super symmetric transformation of the vortex for winding number one unlike those other 
models. 

One way to check the saturation of the bound at the quantum level is to calculate directly 
the quantum corrections in Eq.( |2.17|) . We have shown that the leading quantum correction 



to the mass vanish. The magnetic flux in the right hand side of Eq. ( p,17| ) may not get any 
quantum correction. We still need the quantum correction to the coupling constants in the 
presence of the vortex background to check the validity of the quantum Bogomolnyi bound. 
This is a quite interesting open problem. For some quantities for the particles in the vacuum 
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sector, there exist some perturbative calculations 
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APPENDIX A: SPINS 

In this appendix, we will describe how to determine the spin of elementary excitations. 
The spin of the fermion coupled with the CS field was considered in Ref. ||27|| . It was found 
that the fermions carry spin ±1/2 and the sign of spin is determined by the sign of the mass 
term in the Lagrangian. The spin of vectors in the case of unbroken Maxwell CS gauge 



theory was considered in Ref. while that of broken CS theory considered in Ref. |28j . In 
the following, we will describe how to determine the spin of gauge field of the Maxwell CS 
gauge theory in the broken phase. 

We start from the following Lagrangian 

C = --F^" + -^ X F^A X - \ti 2 A^ (Al) 

This may be considered coming from the spontaneously broken theory of the Maxwell- 
Chern-Simons Higgs in the unitary gauge. The canonical variables of this system are A4 
and 7Tj = F 0i + fe^A 7 . We separate the longitudinal and transverse components using the 
following identification. 

A* = - V iX (A2) 

TTi = eijVjTTp - Vi7T x (A3) 



with the abbreviation V« = Vj/V — V 2 . 
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With these new degrees of freedom, tp, x, ^<p, and ir x , the Hamiltonian density for the 
system in Eq. (|Al| ) is written as 



n 



>7T,, 



K \2 , 1 



xY 



K 2 1 



1 1 



K 



In the derivation of 



, we solved the Gauss law for A as 



By varying the Hamiltonian, we can write the equation of motion as 



a 2 



( v \ 



2 k/_ 
2 



2 K 











) 


[ 7r x/ K J 



with 



(A4) 



(A5) 



(A6) 



K 

2' 



(A7) 



Note that usual notions of variable and conjugate momentum for \ an d tt x has been reversed. 
The mass matrix in the right hand side of ( |A6|) has two eigenvalues. 

-2 / ~2 



2 2 i > 

m ± = + — ± + 



(A8) 



Now, we consider the generaters of the Poincare algebra. The Hamiltonian is the integration 
of its density (|A4|) . The momentum is 



P l = J d 2 x 
The angular momentum is 



J 



1 /Sj /Sj 



(A9) 



(A10) 



This means that ip and 7r x behaves like spin-0 fields. However the boost generators have 
infrared problem which was expected from the decomposition and removing of this infrared 
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divergence will fix the spin, as first noted in Ref. [[| After some manipulation, we have the 
following expression for the boost generator 



B' 



d 2 x 



V 7 - 

x Ti + mpeij—^ip - ipeij 



(All) 



The equation of motion (|A6|) can be rewritten as 



O 2 



I 



9 



\ ( 



V 2 



V 



/ 



\ 



7T-v 



f^)//i 



(A12) 



We can diagonalize the matrix in the right hand side by the new fields £ and rj 



\ 



/ 



(A13) 



///c/iV + (m 2 - Ai 2 )/A^+ 
i (m£ -[i 2 )/N_ iik/N- 
with two suitable normalization constants AT+ and AL. With the new degrees of freedom £ 
and rj the generators of the Poincare algebra have the following form 



H = - I (V.v 



V 



i i 
m_7] 



Pi 
J 

B l 



d 2 x 



£Vi£ + rjVif) 



d OCX Cj 



I 



d 2 x 



V,- V,- 



-V- 



V 



(A14) 
(A15) 
(A16) 
(A17) 



with m± = \Jrn± = y/i 2 + t f' Note that in all formulars, £ and r\ have the same 
contribution except for Bi where the signs of infrared singular terms are opposite. By direct 
application of the argument in Ref. |2| , one can show that after removing the infrared singu- 
larity the angular momentum generator ( |A16| ) has additional terms which are propotional 



to 



\m+\ 

and — 1. 



and 



|m_| 



respectively. From those terms, we can determine spin of £ and rj as +1 



APPENDIX B: SUPERMULTIPLETS 



In this appendix we will describe the supermultiplet structures. We will consider both 
the particle and vortex supermultiplets. The structure will depend on whether we have the 
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central charge or not. In the case of the vortices or the symmetry unbroken sector of the 
vacuum, we will have the central charges while in the particle spectrum in the symmetry 
broken sector we have no central charges. For simplicity we start from the SUSY alsgebra 
without the central charge 

{Q% Ob) = 2af5 AB P' 1 (Bl) 

with 

= <t 3 , *i) (B2) 

We work in the Majorana represention and Q a corresponds to the real spinor in SL(2,R) 
in this appendix. The index A, B denotes the indices for the extended algebra. In the case 
of N — 1, we have only one supercharge Q a . In the massive case (this is sufficient in our 
analysis), we can take the rest frame and the algebra becomes 

{Q a , Q 13 } = 25 al3 M (B3) 

Introducing a = T^giQ 1 - iQ 2 ), ^ = ^jiQ 1 + iQ 2 ) , we get the algebra 

{a, a f } = 1 {a, a} = {a\ a)} = (B4) 

We can construct the Clifford vacuum > with spin j by setting a\Q(j) >= 0. The 

angular momentum operator will then be realized as J = ja^a in this representation. We 
will then have two states, i.e., > with spin j and a^\fl(j) > with spin j + ~. To verify 

this last fact, we note 

[J,Q a ] = \Qa (B5) 

and hence 

[J,a] = ~a, [J,at] = i a t. (B6) 

We now turn to the case of N = 2 without the central charge. The indices A and B will 
run from 1 to 2. In this case, we can rewrite the algebra in Eq . ( [Bl|) as 
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{Q a ,Q* l3 }=2a fl P» (B7) 
with the complex notation Q a = -^(Qi — iQ%)- By defining 

a ^^M iQl + tQ ^ ^^"-^ (B8) 
we get the following algebra; 

{ ai ,a\} = {a 2 ,a\} = 1 (B9) 

{di, a 2 } = {a\, 4} = {ai, at,} = {a 2 , a{} = (BIO) 

and 



[J, ai] = [^, a 2] = ^a 2 . (Bll) 

Based on the algebra, we have four states 

>,a\\n(j) >,al\Q(j) >,a\al\Q(j) >} (B12) 

with angular momentum + |, j — |, j). 

Now consider the supermultiplets in the presence of the central charges. The central 
charge changes the superalgebra 

{Q°, Q*/ 3 } = 2afP» - (a 2 ) Q/3 2C (B13) 

with C = ev $ in our case. In general M > C, and the size of the supermultiplets are the same 
as the massive case without the central charge unless the equlity in the above is saturated, 
where the size of the supermultiplets will be reduced. In the model we are considering, there 
will be the central charges for both topological and nontopological vortices and also for the 
particles in the symmetry unbroken vacuum sector. In all these cases the central charge will 
be equal to the mass. To see how the supermultiplets are reduced in these cases, we define 
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1 



a 2 



2VM 
1 



(Q 1 - iQ< 



1 



i, 



(Q^+tQ 2 *) 



2VM 

(Q 1 + iQ 2 ), at = -^=(Q* 1 - iQ 2 *) 



2VM 2VM 
The superalgebra in Eq.( |B13| ) will then become 

{ai,4} = 0, {a 2 ,4} = l. 



(B14) 



(B15) 



Other anticommutators vanish. The operators ai(a|) should then be realized to be zero. 
Therefore we have only one pair of creation and annihilation operators and hence two states 
(not four) > and a 2 \Q(j) > with angular momentum j and j — |. 

APPENDIX C: PHASE SHIFT ANALYSIS 

In this appendix we will describe the phase shift analysis used for the mass correction. 
First consider k = case. Rewrite the Dirac equation in (|3.10|) as 
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(C2) 



a+j-l/2 _f 
r J 

-f d r + ^ 

We have set /' = \^2evf and dropped the prime for simplicity. From these first order 
differential equations, we get the following 2nd order equations. 



hi(r) 






fh(r) 
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= co 2 






K M r ) 
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\ M r ) 
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(C3) 



and 



Li) 



1 h 2 (r)\ 



h 3 (r) 



I Mr) ^ 

M r ) 



(C4) 



where 
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LtL 



-d 



2 d r +a' 
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-dl 



9r + (j-1/2) 2 



+ / 2 



(C5) 



LLt 



2 _ a-a' , (g+j-1/2) 2 



-dl- 
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(C6) 



Note that LL 1 " can be obtained by changing the sign of a and j in L^L. The equations for 
hi and /i 4 are decoupled. 



-dl- 



d r + a' (a + j + 1/2) 5 



+ 



+ f \h 1 = u 2 hi 



(C7) 



and 



-9! 



9r , (j-1/2) 5 



+ 



f2 ^ - W 2 /l4 



/ 2 )h 4 



(C8) 



The functions will behave near the origin as 



In the asymptotic region of r — > oo, the shape of solutions are 



(C9) 



hi(r) ~ Q!iJ|j- + i/2|(wr) +/?i%+i/ 2 |(wr) 



(CIO) 



and 



/i 4 (r) ~ a 4 J| j _i/ 2 |(wr) +/3 4 iV| j _i/2|(a;r) 



(Gil) 



with Q = \/u 2 — 1. The coefficients of a and /3 will be determined to match the behavior 
near the origin (|C9|). 

The influence of the background vortex will show up as a phase shift of the functions by 
comparing with the case with no vortex background (a = 0,f = 1). The phase shifts are 
given by 



tan5 ? (1) = ^ 

J Oil 



(C12) 
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tan*< 4 > = ^ 



The other two functions h 2 and h 3 are determined by 



( u 


[r) ) 




( u 

hi 
















V) j 




hi 





In the limit r — > oo, we have 



ft 2 (r) 



MO 



9 r /ii + /i 4 
<9 r /i 4 + hi 



5f) + cos(ur - \j - l/2|f 
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First, consider the case with j > 0. The phase shift is given by 

-Q sm(R - 5 { p) + cos(i? - 5f + f ) ^ 
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with 



R = Qr-U + l/2)\-\, K = QT-{j-V2)\-\ = R 



Then the phase shifts become 



tan(5f ) 
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For the case with negative j, the phase shift will be 
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with 
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R" = ur+ (j + l/2)\ - * R"> = ur+ (j - 1/2)| - \ = R" + \ 



The phase shift will then be 



(2) 



tan(5 j - 5 



and 



tan(5f-^ 4) ; 



sinjSf - 5f) 
u + cos(^ 4) - 5f) 
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From the above equations, not depending on the sign of j 



hence 
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From the relation between the phase shift and the density of states EM we get 



n_i_ = n_. 



We turn to the case with k ^ 0. The equations for h,2 and /13 are 
d r + a' (a + j -1/2) 2 
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By solving these two coupled equations, we can get h 2 and h$. hi and are fixed by 

h = - {{d T - -(a + j - \))h 2 - fh 3 ) (C30) 
too V r 2 J 

K = ~/ ~ r ((d r + -(J + h)h - fh) (C31) 
i(uj — k) V r 2 / 

unless a; = and u — k. 

Following the steps for the case of k = 0, we get the result 

5^ + S Xl = 5^ + 5x h (C32) 

which is similar to Eq.( |C27|) 

There may exists a solution to the Dirac equation in (|3.10|) for threshold value of u = k. 
But this gives us a discrete modes which does not contribute to the mass correction. 
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